To APPEAR IN The Astrophysical Journal Letters. 
Preprint typeset using L^TgK style emulateapj v. 26/01/00 



f 



RAPID INFRARED VARIABILITY OF THREE RADIO-LOUD NARROW-LINE SEYFERT 1 GALAXIES: A 
VIEW FROM THE WIDE-FIELD INFRARED SURVEY EXPLORER 

NiNG Jiang'-^, Hong-Yan Zhou' - \ Luis C. Ho'*, Weimin Yuan', Ting-Gui Wang'"^, 
Xiao-Bo Dong'-^, Peng Jiang'-^, Tuo Ji' and Qiguo Tian"* 

To appear in The Astrophysical Journal Letters. 
ABSTRACT 

Using newly released data from the Wide-field Infrared Survey Explorer, we report the discovery of rapid 
infrared variability in three radio-loud narrow-line Seyfert 1 galaxies (NLSls) selected from the 23 sources in the 
sample of Yuan et al. (2008). J0849H-5108 and J0948H-0022 clearly show intraday vai-iabiHty, while J1505H-0326 
has a longer measurable time scale within 180 days. Their variability amplitudes, corrected for measurement 
errors, are ^ 0.1-0.2 mag. The detection of intraday variability restricts the size of the infrared-emitting region 
to ~ 10"^ pc, significantly smaller than the scale of the torus but consistent with the base of a jet. The three variable 
sources are exceptionally radio-loud, have the highest radio brightness temperature among the whole sample, and 
all show detected 7-ray emission in FermifLP^ observations. Their spectral energy distributions resemble those of 
low-energy-peaked blazars, with a synchrotron peak around infrared wavelengths. This result strongly confirms 
the view that at least some radio-loud NLSls are blazars with a relativistic jet close to our line of sight. The 
beamed synchrotron emission from the jet contributes significantly to and probably dominates the spectra in the 
infrared and even optical bands. 

Subject headings: galaxies: individual (SDSS J084957.98H-5 10829.0, SDSS J094857.32H-002225.5, SDSS 
J150506.48H-032630.8) — galaxies: active — galaxies: jets — infrared: galaxies 



quasars 13%). Very radio-loud NLSls, those with radio- 
loudness parameter R > 100, where R is commonly defined as 
the ratio of the radio flux density at 6 cm to optical flux den- 
sity at 4400 A (Kellermann et al. 1989), are extremely rare, 
comprising only 2.5% of the NLSl population (Komossa et al. 
2006). Of particular interest, some individual objects are found 
to exhibit blazar-like behavior (e.g., Zhou et al. 2003, 2005, 
2006; Gallo et al. 2006). Yuan et al. (2008) presented a com- 
prehensive study of 23 genuine radio-loud NLSls (hereinafter 
the Y08 sample). The radio sources in radio-loud NLSls are 
ubiquitously compact, unresolved on scales of several arcsec- 
onds. Some of them show interesting properties, including flat 
radio spectra, compact VLBI cores, very high radio brightness 
temperatures (Tb), enhanced optical continuum emission, flat 
X-ray spectra, large-amplitude X-ray flux and spectral variabil- 
ity, and blazar-like spectral energy distributions (SEDs). Based 
on recent observations taken by Fermi Gamma-ray Space Tele- 
scope (Fermi), some of the radio-loud NLSls indeed display a 
hard X-ray component and prominent 7-ray radiation, which 
undoubtedly arise from a relativistic jets pointed toward us 
(Abdo et al. 2009a, 2009b; Foschini et al. 2011). The blazar 
interpretation of radio-loud NLSls appears to be successful, 
but additional evidence is still needed for the majority of the 
class. Flat and inverted radio spectra are also often seen in some 
normal Seyfert galaxies, although most have much less power 
compared to the objects studied in the YOB sample (e.g.. Ho & 
Ulvestad 2001). 

As a general feature of blazars, intraday variability is ex- 
pected in radio-loud NLSls if the above interpretation is true. 
This has been confirmed by some optical intra-night 
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1. INTRODUCTION 

Active galactic nuclei (AGNs; including Seyfert galaxies and 
quasars), powered by accretion onto supermassive black holes, 
show multiwavelength variability on time scales from years to 
less than a day. As one of the first recognized properties of 
quasars, variability has served as an important tool to investi- 
gate the emission processes of AGNs (e.g., UMch et al. 1997; 
Peterson 2001). In particular, intraday variability can provide 
essential constraints on the central regions of AGNs on phys- 
ical scales smaller than the Solar System, which in general 
cannot be resolved directly by current observational capabili- 
ties. Intraday variability has been detected throughout nearly 
the entire observable electromagnetic spectrum in many radio- 
loud AGNs with flat radio spectra, which are often collectively 
called blazars (see review in Wagner & Witzel 1995). Radio- 
loud AGNs plausibly contain a relativistic jet originating near 
the central black hole. Blazars are believed to be those with 
a relativistic jet oriented close to the line of sight; thus, their 
nonthermal jet emission is highly Doppler boosted. 

Narrow-line Seyfert 1 galaxies (NLSls), originally defined 
by Osterbrock & Pogge (1985), are a special class of AGNs in 
which the broad permitted lines are relatively narrow and the 
[O III] A5007 emission line is weak (see Pogge 2000 for a re- 
view). Apart from this, NLSls also show strong Fe II emission 
lines in their optical/ultraviolet spectra, steep soft X-ray spectra, 
and rapid X-ray variability (Boiler et al. 1996; Leighly 1999; 
Sulentic et al. 2000; Zhou et al. 2006). Zhou & Wang (2002) 
found that the fraction of radio-loud objects among NLSls 
(^ 6%) is significantly less than that found in normal type 1 
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Table 1. WISE Data 



Name 


z 


Wl 


W2 


W3 


W4 




Aw2 


Aw3 


var_flg 


log Ri.4 


logTB 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


J0849+5108 


0.583 


14.03±0.03 


13.14±0.03 


10.04±0.05 


7.54±0.11 


0.11±0.02 


0.19±0.03 


0.15 


4610 


3.16 


13.0 


J0948+0022 


0.584 


13.22±0.03 


12.12±0.02 


9.15±0.04 


6.78±0.08 


0.19±0.01 


0.18±0.02 


0.16 


9903 


2.55 


12.7 


J1505+0326 


0.408 


14.04±0.03 


13.10±0.03 


9.89±0.04 


7.12±0.07 


0.08±0.02 


0.13±0.03 


0.12 


6730 


3.19 


14.0 


others (median) 


0.523 


14.36±0.03 


13.43±0.03 


10.52±0.06 


8.17±0.11 


0.02±0.03 


0.05±0.05 


0.00 


0100 


2.53 





Note. — Col.(l): abbreviated name of object used in this paper and median value of the other 20 sources. Col. (2): redshift. 
Col.(3)-(6): W1-W4 magnitude from WISE All-Sky Source Catalog. Col.(7)-(9): A of W1-W3 magnitudes, where the error is given 
by the median A of standard stars with similar magnitudes. Col. (10): variability flag. Col. (11): logarithm of radio-loudness parameter 
defined as R1.4 = /,/(1.4 GHz)//^ (4400 A), where -R1.4 = 1.9-R for a radio index a = —0.5. Col. (12): logarithm of radio brightness 
temperature (K). Some other related information can be found in Y08. 



ing likelihood of variability (see Hoffman et al. 2012 for de- 
tails). According to varjlg in the WISE catalog, there are tree 
radio-loud NLSls in the Y08 sample with probable variabil- 
ity: SDSS J084957.98H-5 10829.0 (hereinafter J0849H-5108), 
SDSS J094857.32H-002225.5 (hereinafter J0948H-0022), and 
SDSS J150506.48H-032630.8 (hereinafter J1505H-0326), whose 
respective variability flags are "4610," "9903," and "6730." For 
the remaining 20 sources, the variability flags are no larger than 
2 in any of the four bands. The data for these sources are pre- 
sented in Table 1 . 

We first examine the photometric errors of the WISE data 
using standard stars chosen from the SDSS Stripe 82 cata- 
log (Ivezic et al. 2007). The single-frame photometric data 
can be downloaded from WISE All-Sky Single Exposure (Lib) 
Source Table via the GATOR query service of the NAS A/IPAC 
Infrared Science Archive^. We adopt the PSF profile-fit photo- 
metric magnitudes^. The magnitude standard deviation E for a 
given star in a band is calculated as 



observations (e.g., Liu et al. 2009; Shi & Shan 2011). Extend- 
ing such an investigation to the infrared (IR) becomes widely 
possible only after the recent mission of Wide-field Infrared 
Survey Explorer (WISE; Wright et al. 2010), which achieves 
a sensitivity in the 12/xm band more than 100 times higher 
than the most comparable previous mission In frared Astronom- 
ical Satellite (IRAS). Indeed, all 23 objects in the Y08 sam- 
ple are detected by WISE but none of them except one (SDSS 
J163323.58H-471859.0, hereinafter J1633h-4718)0 was detected 
by IRAS. WISE has mapped the entire sky in four bands cen- 
tered at 3.4, 4.6, 12, and 22 ^m (hereinafter the Wl, W2, W3, 
and W4 bands) with angular resolutions of 6"1, 6."4, 6."5, and 
12", respectively. The field-of-view of WISE is 47' x 47' and 
has a small (10%) overlap between adjacent fields in one or- 
bit. The scan circle advances by about 4' per orbit. Thus, 
there are typically 12 successive orbits covering a given source. 
With ^ 15 orbits per day, the observing cadence of WISE is 
well suited for studying intraday variability. It is worth not- 
ing that the perspective from IR variability has some unique 
advantages. First, light in the IR is much less affected by 
dust extinction and/or gas absorption than in the optical, UV, 
or X-ray bands. Even compared with radio, IR is free from 
contamination from refractive interstellar scintillation, which is 
wavelength-dependent in the radio but negligible in the IR. 

This Letter reports the discovery of rapid IR variability in 
three radio-loud NLSls selected from the Y08 sample using re- 
cently released data from WISE. We assume a cosmology with 
Ho = 70 km s-i Mpc"', n„, = 0.3, and SIa = 0.7. 

2. DATA ANALYSIS AND RESULTS 

The WISE All-Sky Data Release Source Catalog contains po- 
sition and four-band photometric data for 563,921,584 objects 
detected on the coadded atlas images (Cutri et al. 2012). Pho- 
tometry was performed by fitting point-spread functions (PSFs) 
simultaneously to all the individual exposures covering an ob- 
ject. A variability flag, varjlg, is assigned to every source 
in each band, giving the probability of flux variation as rep- 
resented by an integer value from to 9. A value of in- 
dicates insufficient or inadequate data to determine the vari- 
ability probability, while values from 1 to 9 indicate increas- 

'However, the IRAS source is most likely associated with the star-forming nuclei of a companion galaxy, separated by 
^http://www.astro. washington.edu/users/ivezic/sdss/catalogs/stripe82. html 
^http://irsa.ipac. caltech.edu/cgi-bin/Gator/nph-dd 

'Brief description of the definition of profile-fit magnitude as well as other related parameters can be found in Cutri et al. (2012). 



(1) 



where (m) is the flux-weighted mean. Figure [T] displays E as 
function of magnitude in Wl and W2 for about 5000 stars. We 
have also overplotted the 23 radio-loud NLSls. The three can- 
didate variable sources from the Y08 sample, as well as another 
one (SDSS J172206.03H-565451.6, hereinafter J1722H-5654), 
have large variance outside the 1-cr scatter of the E-magnitude 
distribution; these outliers are either intrinsically variable or 
their observations have large systematic errors. 

We begin by extracting the light curves to reexamine the vari- 
ability of the three candidate variable objects. The signal-to- 
noise ratio (S/N) is high enough to give reliable photometric 
measurement with S/N > 10 in Wl and W2 and S/N > 5 in 
W3 for all exposures. Judging by the reduced chi-square (x^) 
of the profile-fit photometry, which is close to 1, the PSF fitting 
is reliable and there is little contamination from nearby objects. 
Unfortunately, the W4 data have much lower S/N and larger 
photometric errors; some are marked as "null," which means 
that it is either not measurable or only an upper limit 

4" from the AGN (Yuan et al. 2010). 
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Fig. 1 . — Standard deviation of single-frame pliotometric measurements as a function of magnitude for tlie Wl (left) and W2 (right) band. Black dots are standard 
stars from SDSS Stripe 82; its median values and l-cr upper boundary are marked with the green solid and dashed lines, respectively. The three candidate variable 
sources in the Y08 sample are plotted with blue diamonds, J1722+5654 as a red plus, and the nonvariable sources in our sample as red squares. 



can be placed on the flux. We exclude the W4 data from fur- 
ther analysis. In light of possible photometric zero point offsets 
in different fields, for every source we choose as "standards" 
in the field 30-50 stars with brightness comparable to that of 
the source. For each standard star, we calculate the differential 
magnitude between a single observation and the flux-weighted 
mean magnitude. Then we take the zero point offset for a given 
field as the mean differential magnitude of all the standard stars 
in the field. The offsets in Wl and W2 are tiny, typically <0.02 
mag with a systematic error of ~ 0.01 mag. Most of the stan- 
dard stars are too faint in W3, and thus we could not determine 
a zero point offset for this band 

The variability amplitude is commonly measured by the vari- 
ance of the observed magnitudes, with the contribution from 
measurement errors subtracted. We adopt a formalism similar 
to that used in Ai et al. (2010; see also Sesar et al. 2007), 
whereby 

\ 0, otherwise. 

Here the measurement error ^ includes both the l-cr profile- 
fit photometric error and the systematic error from zero point 
offsets (^zero for Wl and W2, which is, in any case, almost neg- 
hgible): 

1 ^ 

^' = ]^E^' + ^-o- (3) 

Figure 2 shows the final light curves for the three sources for 
which variability has been detected. In summary: 

• J0849+5108 — There are 15 exposures within 1 .2 days, 
yielding Awi = 0.1 1 , Aw2 = 0.19, and Aw3=0.15. Asa 
comparison, the standard stars with similar magnitudes 
in Wl and W2 give a median value of Awi = 0.02 and 
Aw2 = 0.03. 

• J0948+0022 — The 10 exposures within one day give 
Awi = 0.19, Aw2 = 0.18, and Aw3 = 0.16. The standard 



stars with similar magnitudes yield a median value of 
Awi = 0.01 and Aw2 = 0.02. 

• J1505+0326 — This target was observed during two 
epochs separated by about 180 days. No evident in- 
traday variability can be seen in either single epoch 
(Awi,w2,w3 ~ 0), but the magnitude jump between the 
two epochs is considerable (see the bottom panel of Fig- 
ure|2]). The flux-weighted mean Wl and W2 magnitudes 
varied by 0.16 and 0.25 between the two epochs, while 
the observations for both epochs as a whole yield Awi = 
0.08, Aw2 = 0.13, and Aw3 = 0.12. The Wl and W2 
magnitudes are almost equal to those of J0849-I-5108, 
and the A of standard stars is also comparable to that 
for J0849H-5108. Given the redshift of 0.408, the time 
scale of variability should be shorter than 128 days in 
the source rest frame. 

To verify the accuracy of the variability flag assigned by 
WISE, we recomputed the variability probability following the 
formalism devised by Hoffman et al. (2012), but using our 
own error estimates. With the exception of J1633-I-4718, which 
shows likely variability according to our test, our results are 
consistent with those of Hoffman et al (2012). The PSF profile- 
fit magnitudes for J1633-I-4718, however, are probably unreli- 
able because this source is an interacting system of two galax- 
ies separated by ^4" as denoted above, which still cannot be 
resolved by WISE. By contrast, J 1722-1-5654, which shows a 
large magnitude deviation (red plus symbol in Figure [T]), has 
A « after correction for systematic measurement errors. 

Lastly, we performed a cross correlation analysis for the three 
variable sources to search for evidence of time lag between the 
light curves of the three bands. No significant lag was found. 

3. DISCUSSION AND CONCLUSIONS 

Using newly released data from WISE, whose well-designed 
observation mode provides a modest number of repeat obser- 
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Fig. 2. — WISE light curves constructed from profile-fit magnitudes; l-cr error bar are plotted. The red points in the Wl and W2 bands represent measurements 
from a typical nearby "standard" star, shifted slightly by a constant for ease of comparison. 



vations of the sky, we discovered significant IR variability in 
three out of 23 radio-loud NLSls selected from the sample of 
Y08. All three are markedly variable in the Wl (3.4 /im) and/or 



W2 (4.6 /im) bands, at greater than 6a confidence. Closer ex- 
amination of their light curves reveals that J0849-I-5108 and 
J0948H-0022 show clear intraday variability, while J1505H-0326 
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Fig. 3. — SEDs of the three vaiiable objects. Data include radio measure- 
ments (black dots), WISE photometry (red squares), 2MASS near-IR magni- 
tudes (black triangles), SDSS ugriz PSF magnitudes corrected for Galactic ex- 
tinction (black stars). X-ray and 7-ray from Abdo et al. 2009b; D' Ammando et 
al. 2012 (blue diamonds). The mean SED of the other 20 nonvariable sources 
is plotted with a green dotted line. The radio, 2MASS, and SDSS data were 
collected from NED. 

varies on a longer time scale of 180 days or less. Quantify- 
ing the variability amplitude using the variance of the observed 
magnitudes (e.g., Ai et al. 2010), after accounting for mea- 
surement errors, we find A « 0.1 -0.2. This represents the 
first definitive detection of IR intraday variability for radio-loud 
NLSls. 

The variability time scale allows us to place limits on the size 
of the IR-emitting region. In the case of J1505H-0326, for which 
only relatively long-term variability can be determined, the size 
of emission region is confined to be < 0.1 pc. On the other 
hand, the detection of intraday variability in J0849H-5108 and 
J0948H-0022 implies that their IR emission region is very com- 
pact, confined to scales of ^ 10"^ pc, which corresponds to only 
hundreds of Schwarzschild radii assuming Mbh ~ 3 x lO^M© 
(Y08). Thermal IR radiation in AGNs is normally believed to 
be arise mainly from a dusty torus outside of the dust subli- 
mation radius (Laor & Draine 1993), which for graphite grains 
is about 0.5 pc for these three sources. If the inferred size of 
the emission region for J0849H-5108 and J0948H-0022 is rep- 
resentative, it implies that the IR radiation comes from a spa- 



tial scale much smaller than that of the torus, perhaps no larger 
than that of the accretion disk. This strongly suggests that the 
IR emission is dominated by an additional nonthermal com- 
ponent, presumably associated with the jet responsible for the 
radio emission. Coincidently, but perhaps not surprisingly, in 
all three objects variable 7-ray emission has been detected re- 
cently with Fenni/LAT (Abdo et al. 2009a, 2009b; Foschini 
2011; D'Ammando et al. 2012). This confirms the blazar-like 
nature of these three radio-loud NLSls and supports the hy- 
pothesis that they do indeed have relativistic jets pointing close 
to our line of sight. In this scenario, the variable IR emission 
derives from synchrotron radiation from the base of the jet. 

We note, in passing, that intraday variability of J0948H-0022 
with similar amplitude as seen in WISE has been reported also 
in the optical (Liu et al. 2010). IR observations are, in princi- 
ple, superior to optical observations for the purposes of detect- 
ing synchrotron radiation from jets. IR radiation is less affected 
by dust extinction and is less contaminated by AGN accretion 
processes that may cause optical variability due to instabilities 
in the accretion flow. 

The fraction of IR-variable sources in the Y08 sample seems 
to be quite low (3/23 or 13%). One immediate, trivial reason 
is probably simply sensitivity. About half of the Y08 sample is 
fainter than the three variable sources, and their variability, if 
present, is likely masked by measurement errors (cf. Figure 1). 
The degree of radio-loudness also seems to matter The major- 
ity of the parent sample have radio-loudness parameters much 
lower than those of J0849H-5108 and J1505H-0326. J0948H-0022 
has a R parameter 0.6 dex lower than the other two variable 
objects, but it is still higher than half of the full sample (see 
Table 1 and Figure [3]). According to Y08, eight sources dis- 
play prominent radio variability between two epochs separated 
by several years. The three IR variables discussed here show 
the largest amplitude of radio variations (^ 40%-75%) and the 
highest radio brightness temperatures, all exceeding the thresh- 
old value, usually estimated to be 7b ~ 10'^ K, above which 
inverse-Compton catastrophe is predicted to occur (Kellermann 
& Pauliny-Toth 1969). Taking the inverse-Compton limit of 
10'^ K as a conservative limit, the minimum Doppler factors 
are estimated to be 1.5-4.7, which are within the range of val- 
ues for relativistic jets inferred for classical radio-loud AGNs 
(Ghisellini et al. 1993). In addition, they show ubiquitously in- 
verted radio spectra and happen to be the flatest three of the 1 1 
sources whose radio spectral slope have been estimated in Y08, 
which is consistent with the blazar nature. 

Selection by IR variability also must favor objects that are 
intrinsically more IR luminous. The overall SED of blazars 
is characterized by two broad bumps due to synchrotron emis- 
sion and inverse-Compton scattering, respectively. "Blazars" 
is a collective term for BL Lac objects and flat-spectrum radio 
quasars (FSRQs), which are distinguished from each other ac- 
cording to the strength of their emission lines in optical spectra: 
BL Lac objects have very weak emission lines while FSRQs 
look like normal quasars. Depending on the peak frequency of 
the synchrotron bump, BL Lac objects can be divided into low- 
energy-peaked sources, whose SEDs reach a local maximum at 
IR-optical wavelengths, and high-energy-peaked ones, which 
peak near the UV or soft X-ray band (Giommi & Padovani 
1994). Similarly, FSRQs can also be divided into their high- 
frequency-peaked and low-frequency-peaked varieties, parallel 
to the classification devised for BL Lac objects (Perlman et al. 
1998). Our three IR-variable radio-loud NLS 1 s have SEDs that. 
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in fact, peak near the IR (Figure O and qualitatively resemble 
low-frequency-peaked FSRQs in terms of their overall energy 
distribution. (We note that the apparent discontinuity between 
the IR and optical data points on the SED may be caused by 
variability, since the IR and optical data were not taken con- 
temporaneously.) This strongly suggests that our IR selection 
preferentially favors systems with IR-peaked SEDs. 

Two tests support the above hypothesis. First, we exam- 
ined the IR variability properties of the large sample of BL 
Lac objects from Nieppola et al. (2006). Judging from the 
var^flg parameter from WISE, 24 out of the 98 low-energy- 
peaked sources can be considered variable {var^lg > 5) in at 
least one WISE band, whereas, strikingly, only 6 of 110 high- 
energy-peaked source qualify as variable. As in our sample, it 
appears that the higher incidence of IR variability is related to 
the higher IR/optical flux contribution from jets. Further evi- 
dence comes from examining the strengths of the optical emis- 
sion lines. The median equivalent width of H(3 for the three 
IR- variable radio-loud NLSls (23 A) is significantly smaller 
than not only the median value of all quasars (^ 70 A; Zhou et 
al. 2006; Shen et al. 2011) but also that of the parent sample 
of radio-loud NLSls (42 A). This is consistent with the notion 
that IR-variable sources have a stronger optical continuum due 



to enhanced contribution from a jet. 

This study demonstrates that the WISE All-Sky Survey af- 
fords an excellent opportunity to study IR variability in AGNs, 
a powerful diagnostic of their radiation mechanism. Short-term 
IR variability, in particular, is sensitive to spatial scales relevant 
to the base of the jet, providing an effective and efficient probe 
that complements radio and high-energy observations. We will 
take full use of the WISE data to study other classes of AGNs 
in the near future. 
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